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Abstract 

We discuss the dynamics of single particle by laying a hypothesis that the Hamilton's principle of 
stationary action is not exact. We then postulate that the deviation of the action with sufficiently 
short time interval from the stationary action is distributed along a sufficiently long trajectory 
according to an exponential law. We show that the dynamics of the ensemble of trajectories satisfies 
the Schrodinger equation with Born interpretation of wave function if the average deviation is given 
by h/2 and if two opposite signs of deviation occur equally probably. The particle thus behaves as 
if it is guided by a virtual wave satisfying the Schrodinger equation. 
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1 



Let us discuss the dynamics of a particle with mass m in one spatial dimension x subjected 
to a potential V{x, t), where t denotes time. To do this, let us define action along a segment 
of path connecting two spacetime points as follows 

S = j pdx- Hdt, (1) 

where p = mv = d^S and H = —dfS are the momentum and energy of the particle, 
respectively. Here v = dx/dt is the velocity of the particle and we assume that H depends 
on X and p. The Hamilton's principle then says that the only admissible trajectory that 
connects the two points is given by the one in which AS" is stationary Mathematically, 
it is then expressed as 

SS, = 0, (2) 



where the variation is done by keeping the end points fixed, and we have denoted the 
stationary action as Sc{x,t). All the other paths are classically forbidden. 

Performing the variation of Eq. ([5]), one obtains the Hamilton equation dpc/dt = —d^Hc 
and Vc = dx/dt = dp^H^ [1], where we have denoted the classical momentum, velocity and 
energy as pc, Vc and He. Given a specific dynamical problem, the dynamical equation is thus 
determined by choosing a specific form of classical energy He as function of pc and x. In 
particular, the Hamilton equation is equivalent to Newton equation if we choose the classical 
Hamiltonian He as follows: 

-0,Se = He = ^ + V=^^^ + V. (3) 

2m 2m 

Viewed as partial differential equation for Sc{x,t), Eq. is also called as Hamilton- Jacobi 
equation. 

To solve the Hamilton- Jacobi equation one needs to choose the initial stationary action 
Sc{x, 0) which implies an initial classical momentum field Pc{x, 0) = d^Sc- A single trajectory 
is then picked up if one fixes the initial position as well. If one considers ensemble of initial 
position with distribution p{x, 0), then the distribution of the position of the particle at any 
time t is obtained by solving the continuity equation 

dtp = -d^ipVc) = -d^ (p-^] . (4) 

Equation (jlj) conserves the local probability fiow. 
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Let us now exercise an assumption that the Hamilton's principle is not valid exactly but 
retain the deterministic character of the dynamics. Further let us assume that the short 
segment trajectory with stationary action is not the only admissible segment of trajectory 
but the most probable one. This assumption tells us that if we look at a trajectory within 
a large time scale and average short time fluctuations around the stationary action, then 
the trajectory looks effectively classical. However, probing the dynamics within sufficiently 
short time scale will reveal that the momentum and energy are fluctuating around the 
corresponding classical values. 

To proceed one needs to know how the deviation from the stationary action is distributed 
along the trajectory. Let us consider a trajectory running from t = to t = T, slice it 
into sufficiently short segments of equal time interval At = T/N where N is the number 
of segments, and take statistics on the distribution of deviation from stationary action 
IAS* — ASc\ along each segment. Let us then postulate that there is a very short universal 
time interval At such that using At to slice a sufficiently long trajectory so that N is 
sufficiently large, the distribution of \AS — ASd converges into exponential law as 



where n{\AS — ASd) is the number of short segments whose deviation from the stationary 
action is \AS — ASd, h is a constant with action dimension. The average deviation 
is thus h/2. Note that Eq. (j5]) should not be interpreted to give the probability of an 
elementary step of a stochastic dynamics. 

Since our basic law takes statistical form then it is impossible to develop dynamical 
causal relation which refers to single event (trajectory) as in classical mechanics. Hence, 
one is forced to instead consider an ensemble of trajectories described by p{x, t) which is 
transported along the deterministic momentum flow p{x,t) = d^S. Now, let us consider p 
at two very close spacetime points {x,t} and {x + Ax,t + At} connected by a short segment 
of trajectory with an action AS". Since the flow is deterministic and since we are only given 
the distribution of deviation from the stationary action [AS* — ASd along the trajectory, 
then p{x + Ax, t + At) has to be proportional to p{x, t) multiplied by the chance that the 
short segment with deviation IAS" — ASd occurs which is given by Eq. ([5]), and further 
multiplied by a term that describes whether the short segment repels or attracts the nearby 




(5) 



3 



trajectories. The last term thus has to take the form exp{—dxvAt), where d^v = d'^S/m 
is the rate of attraction or repulsion of the nearby trajectories when the sign is negative or 
positive respectively. One therefore has 

p{x + Ax,t + At) ~ p(x, t)e-tl^^-^^=l-^^*. (6) 

Note that Eq. ([6]) does not differentiate between two different cases of AS > ASc and 
AS < ASc. 

Next, assuming that Ax and At are sufficiently small, expanding the exponential on the 
right hand side of Eq. ([6]) up to the first order one has 
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h 
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Ap{x,t) = - - AS- AS, + -diSAt p{x,t). (7) 



m 



Further, expanding Ap and AS* as A/ = ^At + ^Ax, and comparing term by term one 
finally obtains 

^^ = ±{dxS,(x,t)-d,S{x,t)), 
z p 

= ±{dtSc{x,t)-d,S{x,t)) - ^dlS. (8) 

Here the "+" and "— " signs correspond to the case when AS* > AS, and AS' < AS,, 
respectively. It is thus clear that in the regime where the terms containing h in Eq. ([8]) are 
negligible, then the momentum and energy are effectively given by their classical mechanics 
values: 

p = d,S^d,S,, H = -dtS ^ -dtS,. (9) 
As expected, formally, classical mechanics (Hamilton's principle) is regained in the limit 

Before proceeding, let us remark that Eq. ([8]) can not be interpreted as causal dynamical 
relation for single event. The left hand side is determined by probability density p{x, t) which 
gives the relative frequency that the particle is at x at time t in infinitely many trials, and 
the right hand side are dynamical quantities which refer to single event. Equation ([S]) is thus 
descriptive rather than explaining causal relation. Further, to verify the above equations, 
one has to run in principle infinitely many trajectories. As explicitly shown in the original 
equation of ([6]), Eq. (|8]) has to be interpreted as equation for p{x,t) given the dynamical 
variable S{x,t) as in classical physics. Moreover, any causal conclusion about the dynamics 



of the particle based on Eq. (|8]) can only be taken through averaging over the ensemble 
described by p{x,t). 

Now let us discuss how the above statistical modification of Hamilton's principle changes 
the dynamical equation of classical mechanics for ensemble of trajectories. Recall that the 
latter is given by a pair of coupled equations and dl]). Moreover, notice that the coupling 
between p{x,t) and Sc{x,t) in those equations are one sided. Namely, while p{x,t) depends 
on Sc{x,t) through Eq. (jlj), Sc{x,t) does not depend on p{x,t). We shall show that when 
the Hamilton's principle is modified statistically as in Eq. ([5]), then the coupling between 
the action and the probability density will be two sided. 

To do this, first, we shall discuss the case of "+" and "— " signs in Eq. separately 
and combine the resulting equations afterward with additional physical assumption. 

Let us consider the case when AS* > ASc so that Eq. ([8]) takes "+" sign. Inserting the 
upper equation of ([8]) into the continuity equation of (jlj) one obtains 

dtP=-:^dlp-dJp^). (10) 

On the other hand, substituting both equations in ([8]) into the Hamilton- Jacobi equation of 
([3]) one gets 

2m 2m R 

where R = ^/p and we have used the following identity 

^(^y = + _^Mp^ (12) 

8m V p / 2m R 4m p ' 
Finally, inserting Eq. ffTOj) into Eq. ffTTj) one obtains 

H^-a,sJ-^^v-^Lm, (13) 

2m 2m R ^ ^ 

Next, let us consider the case when AS < ASc so that Eq. ([8]) takes "— " sign. Again, 
inserting the upper equation of (IE]) into the classical continuity equation of (jl]) one gets 

d^p=^dlp-dJp^). (14) 
2m V m / 
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Further, substituting both equations in ([8]) into the Hamilton- Jacobi equation of ([3]) one 

gets 

2m 2m R 

where we have used again Eq. ( fT2l) . Substituting Eq. IHM into Eq. ( |T5l) . one finally obtains 

2m 2m R 

which is exactly equal to Eq. flT3|) . 

We have thus two pairs of equations, one is given by Eqs. (fTOj) and (fT3|) if AS* > ASc, 
and another one is given by Eqs. (HM and ( |T3l) if AS < AS'c. In both cases, the energy, 
which is equal to the temporal change of the action H = —dtS, is given by Eq. (|T3l) : whereas 
the temporal change of probability density, dtp, differs only on the sign of the first term on 
the right hand side. Moreover, notice that so far we only assume the statistical distribution 
of the deviations [AS* — ASd of the action. Hence there is still ambiguity in the choice of 
the distribution of cases of AS > ASc and AS < ASc, namely the distribution of the "±" 
signs in Eq. ([H]). Now let us proceed to assume that the two cases of dynamics occur equally 
probably independent of the value of the deviation. Hence, the relative frequency that 
AS > ASc is equal to the relative frequency that AS* < AS'c, namely V{+) = V{—) = 1/2 
regardless the value of [AS* — A^d. Averaging over this fluctuation, the flrst term on the 
right hand side of Eqs. ffTOj) and f|T^ cancel to each other to give 

dtp=-dJp^). (16) 

Finally Eqs. ( |T3l) and (|T6i) can be recast into a compact equation for complex-valued 
function ,..ypexp(,S/fi) as y 

ihdt^ = dlij + Vi). (17) 

2m 

The above equation is just the Schrodinger equation if we identify h = h/{2Tx) where h is 
Planck constant. The Horn's rule is evident \'ip\'^ = p. Further, the ensemble average of 
energy at any moment is given by 

= y dxr[--^di + v'j^. (18) 



The last line is just the quantum mechanical average energy which is conserved by the 
Schrodinger equation of f|T7|) if V is independent of time. 

To conclude, we have discussed the dynamics of ensemble of trajectories of single particle 
by assuming that the Hamilton's principle of stationary action is not exact. We then pos- 
tulate that there is a universal short time interval so that the deviation of the action from 
the stationary action is distributed along the trajectory according to exponential law with 
average h/2. The dynamics of the ensemble of trajectories is then shown to be governed 
by the Schrodinger equation with Born interpretation of wave function if the two opposite 
signs of the deviation occur equally probably. 

Let us compare the dynamics of ensemble of trajectories developed in this paper with 
Nelson stochastic dynamics ^ and de Broglie-Bohm pilot- wave theory two dynamical 
theories for ensemble of trajectories which are also governed by the Schrodinger equation. In 
stochastic dynamics, the particle is assumed to interact with some universal background field 
so that it undergoes a stochastic Brownian-like trajectory with diffusion constant h/{2m). 
Hence the dynamics is stochastic. In contrast to this, the dynamics in our model is determin- 
istic. Moreover, in stochastic dynamics, it is the trajectory which is subject to fluctuations, 
whereas in our dynamical model, it is the action itself which fluctuates. 

On the other hand, similar to our dynamical model, the pilot-wave theory is determinis- 
tic. However, in contrast to pilot-wave theory in which the wave function is assumed to be 
physically real, in our dynamical model, the wave function is merely an artificial mathemat- 
ical tool to describe the ensemble of trajectories of a single particle. Hence, borrowing the 
language of pilot-wave theory, in our model, the particle behaves as if it is guided by a vir- 
tual wave, tpi^jt), satisfying the Schrodinger equation of (fT7|) . This leads us to expect that 
similar to the pilot-wave theory, our model of dynamics will also show statistical wave-like 
pattern in double slit experiment {g]. 
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